arXiv: 1502.03693v3 [cond-mat.soft] 14 Sep 2015 


Dynamic scaling for the growth of non-equilibrium 
fluctuations during thermophoretic diffusion in 
microgravity 

Roberto Cerbino^ Yifei Sun^, Aleksandar Donev^, and Alberto Vailati^ '^ 

'Universita degli Studi di Milano, Dipartimento di Biotecnologie Mediche e Medicina Traslazionale, Milan, 1-20133, 
Italy 

^New York University, Courant Institute of Mathematical Sciences, New York, NY 10012, USA 
^Universita degli Studi di Milano, Dipartimento di Fisica, Milan, 1-20133, Italy 
*albe rto.vailati@unimi.it 

ABSTRACT 

Diffusion processes are widespread in biological and chemical systems, where they play a fundamental role in the exchange 
of substances at the cellular level and in determining the rate of chemical reactions. Recently, the classical picture that 
portrays diffusion as random uncorrelated motion of molecules has been revised, when it was shown that giant non-equilibrium 
fluctuations develop during diffusion processes. Under microgravity conditions and at steady-state, non-equilibrium fluctuations 
exhibit scale invariance and their size is only limited by the boundaries of the system. In this work, we investigate the onset 
of non-equilibrium concentration fluctuations induced by thermophoretic diffusion in microgravity, a regime not accessible 
to analytical calculations but of great relevance for the understanding of several natural and technological processes. A 
combination of state of the art simulations and experiments allows us to attain a fully quantitative description of the development 
of fluctuations during transient diffusion in microgravity. Both experiments and simulations show that during the onset the 
fluctuations exhibit scale invariance at large wave vectors. In a broader range of wave vectors simulations predict a spinodal-like 
growth of fluctuations, where the amplitude and length-scale of the dominant mode are determined by the thickness of the 
diffuse layer. 


Introduction 

Diffusion in liquid mixtures and suspensions represents a fundamental spontaneous mass transfer mechanism at the microscopic 
scale. For instance, it regulates transport processes in the cell, the growth of crystals and the kinetics of aggregation 
of macromolecules and colloidal particles in suspension. During the last 20 years it has been shown, both theoretically 
and experimentally, ' that diffusion is accompanied by non-equilibrium concentration fluctuations exhibiting generic scale 
invariance- in the length scale range from the molecular scale up to the macroscopic size of the system. Quite interestingly, it 
has been shown that non-equilibrium fluctuations do not represent merely a perturbation of a macroscopic state; instead, the 
diffusive flux can be understood to be entirely generated by non-equilibrium fluctuations. On Earth the scale invariance of 
the fluctuations is broken at small wave vectors by the presence of the gravity force^ that either quenches^’ ^ or amplifles‘^“* ' 
long wavelength fluctuations, depending on whether the density profile is stabilizing or not. Under microgravity conditions the 
scale invariance is broken by the finite size of the diffusing system.*^ 

So far, theoretical models suitable to describe the statistical properties of non-equilibrium fluctuations have been developed 
only for systems at steady state' or for systems whose macroscopic state evolves much slower than the fluctuations.*^ While 
this is always the case for weakly confined systems undergoing diffusion on Earth, the situation in microgravity conditions is 
more complex and no theoretical model is currently able to provide an analytic description of non-equilibrium fluctuations 
occurring during transient diffusion processes in microgravity. This is due to the fact that under microgravity conditions 
the modes associated to the macroscopic state and to the fluctuations evolve with the same timescales, thus preventing the 
separation of the two contributions. However, since the advent of space platforms, several experiments controlled by diffusion 
have been performed in a microgravity environment, which guarantees the absence of spurious convective motions. Notable 
examples include experiments on the crystallization of proteins,*"*’ *^ critical phenomena,*^’ *' the investigation of the influence 
of vibration on diffusion,*^’ *^ and of transport properties in ternary mixtures.-** Therefore, the understanding of the onset of 
concentration fluctuations during diffusion in the absence of gravity represents an important feat both from the fundamental 
point of view, due to the lack of suitable theoretical models, and from the experimental point of view, due to the huge investment 
required to perform experiments in Space. 


In this work, we investigate both experimentally and computationally the onset of non-equilibrium concentration fluctuations 
in a polymer suspension under microgravity conditions. We quickly apply a temperature gradient to the initially homogeneous 
polymer solution. The gradient gradually induces the formation of a concentration profile through thermophoresis (Fig. 1).“’ 
The time evolution of the fluctuations is monitored experimentally by using a quantitative shadowgraph technique.^“’^^ The 
fluctuations are also simulated under the same conditions by using a finite-volume method recently developed for the study 
of giant fluctuations in confinement.^'*’"'’ For large wave vectors, the scale invariance of the fluctuations is confirmed, both 
by experiments and simulations, also during the transient. Interestingly, simulations predict that a dominant mode in the 
structure factor of the fluctuations is found at small wave vectors during transient diffusion. The wave vector associated to 
this dominant mode decreases as time goes by, with a kinetics compatible with a diffusive growth. For long time, the peak 
disappears and is replaced by the expected plateau due to the effect of the impermeable boundaries.The kinetics observed 
during the transient bears many similarities with that of spinodal decomposition,"^ the most notable feature being that the 
structure factor S{k,t) of the fluctuations at different times t can be scaled onto a single master curve F{k/km) by using a 
scaling relation S{k/km,t) — km{t)^°‘F{k/km)^^~^^ 
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Figure 1 . Numerically calculated time evolution of the concentration profile 


Results 

Experiments have been performed aboard the FOTON M3 spaceship by using the GRADFLEX facility developed by ESA. 
Eoton M3 is an unmanned spaceship orbiting at an average distance from the Earth of the order of 300 km. The great advantage 
of such a platform with respect to other facilities, such as the International Space Station, is the very low level of residual gravity, 
of the order of 0.7 /rg on average. The GRADELEX setup comprises a thermal gradient cell and a quantitative shadowgraph 
optical diagnostics (Fig. 2). The sample is a suspension of polystyrene (molecular weight 9100) in toluene with a weight 
fraction concentration of 1.8%. It is contained inside the thermal gradient cell’"’ *^ whose thermal plates are two sapphire 
windows. These windows are in thermal contact with two annular thermo-electric devices. This peculiar configuration of the 
cell enables using the sapphire windows both as thermal plates and as observation windows for the detection of fluctuations. 
The temperature of the sapphire windows is monitored by using thermistors that drive two Proportional-Integral-Derivative 
servo control loops. The light source is a superluminous Light Emitting Diode coupled to an optical fiber. The collimated light 
coming from the diode crosses the sample, where it gets partially scattered by non-equilibrium fluctuations. The superposition 
of the scattered light and of the main beam gives rise to an interference pattern onto the sensor of a Charged Coupled Device 
camera. This pattern can be analyzed statistically by using the theory of quantitative shadowgraphy-"’^* to determine the 
structure factors of temperature and concentration non-equilibrium fluctuations. A typical measurement run performed in 
space involves the automated execution of a stabilization phase followed by a measurement phase: after a stabilization of the 
equipment lasting 3 hours, the sample is kept at a uniform temperature T — 30"C for 90 minutes; the quick imposition of a 
temperature difference AT at time f = 0 determines the start of the diffusive process. The typical time constant associated to the 
growth of temperature difference across the sample is of the order of Xj « 100s, significantly smaller than the time needed for 
the diffusive process to reach a steady state Xc ~ 2000s. The presence of a temperature gradient inside the sample determines 
a non-equilibrium thermophoretic contribution to the average mass flux j = —pD[Vc-|-57’c(l — c)V7’] that gives rise to the 
development of an almost exponential concentration profile (Fig. 1).^" Here D = 1.97 x 10^®cm^/s is the mass diffusion 
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Figure 2. GRADFLEX setup. The sample is sandwiched between two sapphire windows in thermal contact with 
thermo-electric devices. The light coming from a Light Emitting Diode through a hber is steered by a mirror and collimated 
onto the sample by a lens. A relay lens collects the light from the main beam and the light scattered by the sample, which 
interfere onto the sensor of a CCD camera. The light path is kept under vacuum. 


coefficient and St = 6.49 x 10 * is the Soret coefficient. At steady state and at the impermeable boundaries the net flux 

must vanish and the concentration profile is characterized by a gradient Vc = —57'c(l — c)V7’. 

Non-equilibrium temperature and concentration fluctuations at steady state 

The simultaneous presence of a temperature and a concentration gradient determines the onset of both temperature and 
concentration fluctuations. For the polymer suspension of interest here the diffusion coefficient is much smaller that the 
thermal diffusivity K — 8.95 x lO^'^cm^/s of the sample, and the timescales for the relaxation of temperature and concentration 
fluctuations are well separated. The wide difference between these timescales was used in Ref. 12 to estimate the Fourier power 
spectrum of concentration fluctuations by using a standard dynamic analysis.Here we use a refined procedure that allows 
obtaining also the power spectrum of the temperature fluctuations. In addition, we use the power spectra of both temperature 
and concentration fluctuations to estimate the corresponding structure factors at steady state (Fig. 3). The advantage of this 
procedure lies in the fact that it allows a precise determination of the temperature difference AT across the sample, which 
was not measured directly in the GRADFLEX experiment but rather estimated from thermal modeling of the sample cell. 
Eitting (bottom dashed line in Eig. 3) the temperature SgQ{k) to the analytical theoretical expression determined by De Zarate 
and Sengers by using a Galerkin approximation’’^^ provides the estimate AT = 13.25K, which is 24% smaller than what was 
previously estimated by thermal modeling.’- The Galerkin approximation systematically under-estimates the structure factor 
at small wave numbers^"^ by a factor of 500.5/720 = 0.695 and is therefore a source of additional error; the computational 
method used here does not make any such uncontrolled approximations. There is presently no exact closed-form theoretical 
expressions available for perfectly conducting boundaries. 

Once a reliable estimate for AT has been obtained, we applied a similar procedure to obtain the structure factor of 
concentration fluctuations at steady state in absolute units. The experimental estimate turns out to be systematically slightly 
smaller than the theoretical predictions made by using a recent exact prediction-^"’ (solid line in Fig. 3). We believe that this 
discrepancy can be attributed to an actual concentration of the sample about 10 % below the nominal value of c = 1.8 % w/w. 
The results of simulations are also shown in Fig. 3 (dashed lines) and at large wave vectors are in fair agreement with both 
experiments and theory. In Fig. 3, it can be noticed that we could not obtain experimental results at very small wave vectors. 
This is due to the presence of a drift of the optical background of the shadowgraph setup for long times, which prevents the 
characterization of the concentration fluctuations at small wave vectors, but in principle does not affect much the short-lived 
temperature flucmations. 

At large wave vectors, the structure factors of both temperature and concentration fluctuations scale as mirroring the 
scale invariance of the fluctuations. However, at a wave vector kf^ k, n/h, the finite thickness h of the sample along the applied 
gradient produces different effects on the two structure factors because of the different boundary conditions for concentration 
and temperature. Indeed, the boundaries are impermeable to mass but conduct heat very well. As a consequence, long wave 
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Figure 3. Structure factors of the non-equilibrium temperature (data on bottom and right y-axis) and concentration (data on 
top and left y-axis) fluctuations at steady state under microgravity conditions. Circles: experimental results; dashed lines: 
simulations; solid lines: theory. 


length temperature fluctuations can be dissipated effectively through the boundaries and a peak in the temperature Sgg{k) 
can be observed. In contrast, in the case of concentration fluctuations the boundaries are impermeable and long wavelength 
fluctuations can be dissipated by diffusion only, which leads to a plateau in Scc{k) for 

Onset of non-equilibrium concentration fluctuations 

The selected experimental sample represents an ideal system to investigate the onset of concentration fluctuations. In fact, 
the small diffusion coefficient determines the progressive development of a macroscopic concentration profile lasting about 
30 minutes. The sample is initially kept at a uniform temperature of 30‘’C. The diffusion process is started by imposing a 
temperature difference AT — 13.25 K at f = 0. Every 10 s we record a shadowgraph image of the sample. The long timescale 
associated to the development of a macroscopic concentration profile enables us to grab 200 shadowgraph images of the sample 
during the approach to steady state. Due to the fact that the system is evolving in time during the transient, it is not possible 
to recover the structure factors of non-equilibrium fluctuations during the transient by applying the same procedure used to 
recover them at steady state. Instead, in this case we rely on a dedicated processing algorithm that takes advantage of the fact 
that after about 100 s the temperature profile reaches a steady state. For this reason, the first 190 s of the process have not been 
included in the analysis. Starting from the image taken at f = 200 s, structure factors of the concentration fluctuations have 
been averaged on groups of 10, 15, 20, 30, 80 images, corresponding to average times of 245, 370, 545, 1345 s. This procedure 
allows reducing the noise on the structure factor by increasing the statistical sample, without losing much temporal resolution. 

Due to the lack of a theoretical model suitable to deal with a transient system, we have performed simulations under 
conditions and sampling procedure mirroring those found in the experiment. A comparison of the experimental and simulated 
data is shown in Fig. 4. The experimental results are in fair agreement with those of the simulations, the only exception being 
the small k behavior of the structure factors corresponding to 245 s and 370 s. For such times, an effective subtraction of the 
optical background is difficult due to the sudden application of the temperature difference, which is particularly limiting when 
the signal is weakest. To partially avoid these disturbances the optical path is kept under vacuum, but when the light scattered 
by the fluctuations is weak the signal at small wave vector is dominated by fluctuations in the optical path of the probe beam 
and by mechanical drifts of the setup. This effect limits our accessible range and prevents the experimental observation of 
a peak in the structure factors, which is well visible in the simulation results only during the short-time transient. This peak 
is associated to the presence of a dominant mode with a wave vector that gradually decreases in time (Fig. 5a), while the 
amplitude of the mode increases progressively (Fig. 5b). 

A first understanding of the presence of a peak can be achieved by taking into consideration that in the presence of fully 
developed temperature and concentration profiles the structure factor S{k,t) of concentration fluctuations grows diffusively. In 
fact, under these conditions it can be shown that, ignoring the influence of the boundaries, S{k,t ) oc [1 — exp{—2Dk^t)]S{k,°o) 
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Figure 4. Time evolution of the structure factor of non-equilibrium concentration fluctuations during the approach to steady 
state under microgravity conditions. Solid lines; experimental data; dashed lines: simulations; dashed-dotted line: exact 
theory’^ 


where S{k,°o) is the structure factor at steady state. A simple model along this lines provides the right qualitative behavior and 
gives rise to a peak in the structure factor behaving asymptotically as and at small and large k, respectively. However, in 
any real diffusive process the modes associated to fluctuations and to the macroscopic state evolve with the same time constants. 
Therefore, the assumption of an initial fully developed concentration gradient on top of which fluctuations develop is rather 
unrealistic. In practice, real effects like the progressive development of a temperature gradient and the subsequent growth 
of boundary layers in the concentration profile are difficult to model theoretically, but can be modeled exactly by means of 
simulations. The fit of the peak of the structure factors of the simulated fluctuations allows us to recover the wave vector k,„ and 
the structure factor Sm = S{km) of the dominant mode. The time evolution of k„, at times smaller than about 100 s is compatible 
with a diffusive growth of the mode k„, (Fig. 5a). During this initial phase, the two boundary layers grow without 

feeling much the presence of each other. However, after a time T = {hjT)^ j(n^D) « 120s they reach a thickness comparable to 
hjl, and the system enters a diffusive regime where finite size effects become relevant, as mirrored by the slowing down of the 
decrease of k„,. 

Dynamic scaling of non-equilibrium concentration fluctuations 

An important feature of the dynamics of the dominant mode is the relation between and Sm (Fig. 5b). One can appreciate 
that at the times larger than 200 s, when the system has entered the restricted diffusion regime, km and Sm are related by a power 
law km Sm, with exponent j3 = 0.12. 

Qualitatively, this behavior is similar to that reported for spinodal decomposition"^"^*' and other phenomena, such as colloidal 
aggregation.^^ The growth dynamics of the structure factor of the concentration perturbations associated to these phenomena is 
such that the structure factor exhibits dynamic scaling S{k/km,t) = k"“F’(k/km), where F{k,km) is a time independent master 
curve and a a power law exponent, which in the case of spinodal decomposition corresponds to the dimensionality of the 
space. This suggests that our results are compatible with a scaling law akin to that of spinodal decomposition with a power 
law exponent a = l/j3 « 8. By scaling the structure factors of simulations in the time range 200i <t< 2QQQs using the 
relationS{k/km,t)k^ = F{k/km) we get that the curves nicely collapse onto a single, time-independent, master curve F{k/km) 
(Fig. 6). 
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Figure 5. a) Time evolution of the wave vector of the dominant mode. The dashed line corresponds to a diffusive behavior 
km b) Wave vector of the dominant mode as a function of the mean squared amplitude (power) of the mode. The 

— 1/8 

dashed line correspond to a power law behavior km Sm 


Discussion 

It turns out that there are some qualitative analogies between the growth of non-equilibrium concentration fluctuations and 
that of the domains in spinodal decomposition. In the case of spinodal decomposition the presence of a dominant mode is due 
to the fact that the process is controlled by a generalized diffusion equation where the diffusion coefficient is negative. This 
uphill diffusion determines the growth of the domains and the progressive buildup of large concentration gradients. In the case 
of non-equilibrium fluctuations we know that in microgravity both the macroscopic state and the fluctuations are controlled 
by a diffusion equation with positive D in the presence of a steady counter-flux determined by the Soret Effect. At steady 
state the diffusive and Soret fluxes balance each other and there’s no net mass flow through the sample. However, during the 
transient the mass flux is dominated by the Soret contribution, and the net balance in the flux of mass determines the growth of a 
concentration gradient, similarly to what happens during the demixing process that drives spinodal decomposition. For spinodal 
decomposition the power law exponent used for the scaling is the dimensionality of the space; in our case it is close to 8. 

Our results provide experimental evidence that linearized fluctuating hydrodynamics quantitatively describes the time- 
dependent growth of fluctuations during transient diffusion processes. Our experimental results are calibrated and compared to 
theory in absolute units, thus significantly extending previous studies for steady-state fluctuations in microgravity. Analytical 
calculations are essentially infeasible in the presence of a transient reference state, especially in the absence of separation 
of time scales as in diffusive mixing in microgravity. The development of numerical techniques for solving the equations of 
fluctuating hydrodynamics"'*’^^ has allowed us to predict the existence of a dynamic scaling law during the development of 
non-equilibrium fluctuations that has not yet been observed in experiments. 


Methods 

Measurement setup 

The GRADFLEX Mixture setup comprises a thermal gradient cell and a shadowgraph optical diagnostics. The gradient cell*"’-^^ 
consists of two - 12 mm thick - sapphire windows kept at a distance of 1.00 mm from each other by means of a calibrated spacer. 
The lateral confinement of the sample is achieved by means of a Viton gasket with an inner diameter of 27 mm. The relatively 
high thermal conductivity of sapphire guarantees a temperature uniform within 3% across the contact surface of the window 
with the sample. Each sapphire window is sandwiched with a coupling ring made of aluminum that brings it into thermal contact 
with an annular thermo-electric device with an inner bore with a diameter of 27 mm. The temperature of the sapphire windows 
is monitored by Negative Temperature Coefficient thermistors that drive two independent Proportional-Integral-Derivative 
servo-controls, which allow to achieve a stability of the temperature of the windows the order 10 mK over 24 hours. The optical 
shadowgraphy diagnostics makes use of a super-luminous Light Emitting Diode with a wavelength of 680 nm and a bandwidth 
of 13 nm. The LED is coupled to a mono-mode optical fiber. The diverging beam coming out of a fiber is steered by a mirror 
and collimated by an achromatic doublet. The role of the steering mirror is to fold the optical path, to maintain the size of 
the instrument compact. The collimated beam goes through the sample, where it gets partially scattered by non-equilibrium 
fluctuations. The main beam and the scattered light are collected by a relay lens and superimposed onto the sensor of a Charged 
Coupled Device camera with a resolution of 1024 x 1024 pixel and a pixel depth of 10 bit, which records an image every 10s . 
In order to avoid disturbances generated by air, the light path is kept under vacuum by means of vacuum tube which can be 
connected to the outer environment of the spaceship by means of a remotely actuated valve. 
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Figure 6. Dynamic scaling of the spectra of the non-equilibrium fluctuations during the approach to steady state. The inset 
shows the unsealed structure factors. Times span the range 200i St = 2000i and are distributed geometrically with a 
multiplier of 1.21 (for a total of 13 curves). 



Optical diagnostics 
Quantitative shadowgraphy 

The non-equilibrium temperature and concentration fluctuations arising as a consequence of the application of a macroscopic 
temperature gradient to a polymer solution (from here-on the sample) give rise to refractive index fluctuations that can be 
detected by using optical shadowgraphy.The phase of a plane wave of intensity Iq that impinges on the sample is 
locally altered by any refractive index inhomogeneity, which causes light scattering. Quantitative shadowgraphy is based 
on the idea that, sufficiently far away from the sample and for weakly scattering systems, the scattered light interferes with 
the transmitted plane wave creating a time-dependent hologram I(x,y,t) ~ /q + 2v^Re[£’i(x,y,f)] at some distance z from 
the sample.The shadowgraph signal is defined as s{x,y,t) = [I{x,y,t)/Io] — 1 = 2Re[Es{x,y,t)/Eo], where Es{x,y,t) is the 
amplitude of the electric field scattered from the sample at distance z- If we indicate with s{k,t) the spatial two-dimensional 
Fourier transform of s{x,t) (from hereon k = {kx,ky) and x = (jc,y)), then the Fourier power spectrum of the shadowgraph 
signal is given by 


A(k) = 



= AVTF{k)kl 


iMf 


— Acc(k) +AgQ(k) 


( 1 ) 


where V is the imaged volume, ko is the wave-vector of the incident light, n{c,T) is the refractive index, 7> is the transfer 
function of shadowgraphy, and Sec and Sgg are the structure factors of concentration and temperature fluctuations as defined in 
Ref. 1 , respectively. Knowledge of the transfer function 7> (k) is thus needed for the quantitative determination of the structure 
factors of the fluctuations and it requires a suitable calibration of the optical setup. 


Transfer function caiibration 

The optical setup was calibrated by using polystyrene spheres with a nominal diameter of 2.0 pm, dispersed in isopropyl 
alcohol.'- Such sample provides a large optical contrast that in the wave-vector range accessible to our experiments gives rise 
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to a constant scattering intensity, representing thereby the ideal calibration sample. The amplitude A„,/(^) determined by using 
the calibration sample was fitted to the function Acai(k) = U(k) sin^ [k^z/{2ko) + 0] +V(k), where k = ^k^ + k^ and where 

(j) = 1.78 was found to match expectations from Mie theory.The so determined U and V were thus used to reconstruct the 
transfer function Tp for the fluctuations that is obtained when 0 = 0 is used. Once the transfer function Tf (A:) is known, the 
setup can be used for the quantitative assessment of the static and dynamic scattering properties of the sample. 


Reduction of experimental results 

Dynamic analysis and isolation of the concentration and temperature contributions at steady state 

A typical analysis of the shadowgraph images I{x,y,t) acquired at steady state at various times t involves the process¬ 
ing of 5000 images. By using a variant of the differential dynamic algorithm the image structure function Dj{k,At) = 

I{k,t + At) —I{k,t) \ is calculated by averaging over pairs of images separated by the same Here l(k,t) is the 


spatial two-dimensional Fourier transform of I{x,t). Theoretical expectation is that Dj{k,At) = 2A(k) 1 — f(k,At) +2B{k), 
where A(k) is given in Eq. 1, B(k) is a dynamic background term that accounts for the noise in the detection chain and 


f(k,At) is the intermediate scattering function of the fluctuations. For our experiments f(k,At) = , 

A(A:) A{k) 

where Tc(k) and Trik) are the characteristic correlation times of the concentration and temperature fluctuations, respec- 

Acc{k)(^-e^^n+Ae9{k)(l-e^^) 


tively. One thus has Di[k,At) = 2 


^2B(k). In practice, in our k-range TT{k) 


is smaller than the time elapsed between the acquisition of two successive images (10 s). As a result, temperature fluc¬ 
tuations appear as uncorrelated background signal and their static scattering contributes to the background. One has 

Dj(k,At) = 2Acc(k) —e + 2Beff(k), where the effective background Beff(k) = B(k) +AgQ(k) incorporates also 

the static scattering signal Agg (k) from temperature fluctuations. Fitting of the experimental curves for Dj(Jc,At) provides thus 
estimates for Acc(k), Beff{k) and Tc(k). We also independently determined B{k) from the dynamic analysis of images acquired 
in the absence of any temperature and concentration gradients, which in turn enabled us obtaining also an estimate for Age (k). 
Using Eq. 1, we recovered the structure factors Scc(k) and Sgg(k) at steady state (Fig. 3). 

Analysis of the transient 

The 200 shadowgraph images acquired during the transient contain contributions coming from concentration fluctuations, 
temperature fluctuations, dynamic noise and static background. The Fourier power spectmm B{k) of the dynamic background is 
time independent and can be characterized accurately from the dynamic analysis at steady state, as described above. Similarly, 
after a time of about 100s needed for the onset of temperature fluctuations, their contribution Agg(k) to the Fourier power 
spectrum of the shadowgraph signal becomes time independent and coincides with that determined at steady state. Conversely, 
the static background contribution represents a time-independent additive term to each shadow image arising from a non-uniform 
illumination of the sample. This contribution could be eliminated easily by using the differential dynamic analysis described 
above. However, the differential dynamic analysis cannot be applied to the images taken during the transient due to the limited 
statistical sample. To overcome this limitation, we determined the static background term by averaging in time 700 images 
collected at steady state IsBix,y) = {I{x,y,t)). The shadowgraph signal during the transient is then defined as Str{x,y,t) = 

[I{x,y,t)/IsBix^y)] —\ and its Fourier power spectrum is given by A,r (k) = ^ Str(k,t) ^ = Acc,tr(k) +Agg(k) + B(k). The 

temporal average was processed by skipping the first 19 images, to avoid effects related to the onset of temperature fluctuations, 
and by averaging the following images in groups of 10, 15, 20, 30, 80, corresponding to average times of 245, 370, 545, 1345 s, 
respectively. The structure factor of transient concentration fluctuations in absolute units was then be determined from the 

relation S(k) = Atrik) — AggCk)-B{k)\/ 




Numerical simulations 

We performed computer simulations of the experimental setup using finite-volume methods for fluctuating hydrodynamics 
described in more detail elsewhere;-^’"'’’^’ here we summarize some key points. In particular. Section V.A of the work of 
Delong et al.-"* present simulations of giant fluctuations in the GRADFLEX experiment, which are used as a basis for the 
more detailed computations reported in this work. The numerical methods have been implemented in the IBAMR software 
framework.'*- We will use CGS units in what follows. In the numerical computations we align the gradient with the y axes in 
order to unify the notation for two and three dimensional simulations. 
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Our numerical codes solve the following stochastic partial differential equations for the fluctuating fluid velocity field 
v(r,t), the mass concentration c(r,f), and the temperature T 

pd,v + Vn =T 7 V^v +V- (2) 

V-v=0 

4c + yVc=DV-(Vc + c(l-c)57’V7’) (3) 

dtT + vVT =kV^T, (4) 


where (r,f) denotes white-noise stochastic forcing driving the thermal fluctuations in the momentum flux (stochastic 
stress). Here 77 = 5.21 • 10^^ is the shear viscosity, n {r,t) is the mechanical pressure. To = 298 is the average temperature, 
St = 6.49 • 10^^ is the Soret coefficient, D=\.91 ■ 10^® is the diffusion coefficient, and tc = 8.95 • 10^^ is the thermal diffusivity. 
The boundary conditions for the velocity are no-slip on the bottom and top sapphire walls, while the other directions are 
periodic. We will discuss boundary conditions for temperature and concentration shortly. 

A number of physical approximations have been made in formulating the system of equations (2,3,4). First, we have ignored 
thermal fluctuations in the mass flux and in the heat flux, which are responsible for equilibrium fluctuations in the concentration 
and temperature; this is justified since our focus is on the much larger non-equilibrium fluctuations. Second, we have used a 
constant temperature Tq for the stochastic stress tensor instead of a spatially-varying temperature; this is justified because the 
maximum difference in temperature across the sample is on the order of a tens of degrees. Third, the density p = 0.858 is taken 
to be constant in a Boussinesq approximation. 

In linearized fluctuating hydrodynamics the equations (2,3,4) are expanded to leading order in the magnitude of the 
fluctuations 5c = c — (c), 5T = T — (T) and 5v = v — (v) = v around the steady state solution of the deterministic equations.' 
As explained in detail in Ref. 24, our numerical methods perform this linearization numerically by solving 

the fully nonlinear equations with weak noise. For the example studied here, in the linearized fluctuating hydrodynamics 
regime, there is no difference between two and three-dimensional simulations due to the symmetries of the problem. Also note 
that in microgravity the temperature and concentration fluctuations are completely decoupled since there is no buoyancy terms 
feeding back into the momentum (velocity) equation. Therefore, numerically we separately solve (2,3) for concentration when 
examining concentration fluctuations, and we separately solve (2,4) when examining temperature fluctuations. The reason 
for this is that these two cases require different temporal integrators, as explained in extensive detail in Ref. 24. We therefore 
separately discuss concentration and temperature fluctuations. 

The experimentally observed light intensity, once corrected for the optical transfer function of the equipment, is proportional 
to the intensity of the fluctuations in the concentration and temperamre averaged along the gradient. The contribution due to 
concentration fluctuations to the shadowgraph is therefore related to the Fourier transform cx {k,t) of the vertically averaged 
concentration, CT{x,z',t) = Jq c{x,y,z',t)dy, where L = 0.1 is the thickness of the sample. More specifically, our simulations 

compute the time-dependent static stmcture factor S (kjc,k^;t) = / f^cx) f^cx) ') , and similarly for temperature fluctuations. 


Concentration fluctuations 

Typical liquid mixtures have a large Schmidt number. Sc = V/D ^ 1, in particular, for the GRADFLEX mixture « 3 • 10^. 
This makes direct numerical solution of the original inertial equations (2,3) numerically infeasible; the time step size needs to 
be chosen to resolve vorticity fluctuations but the time scale of interest is the much longer mass diffusion time scale. Therefore, 
we first take a limit of equations (2,3) as Sc —t in the linearized setting this overdamped limit amounts to deleting the inertial 
term pd/V in the velocity equation.Lastly, it is convenient to approximate the Soret flux c (1 — c) with the linearization 
cSt, which is valid since c ^ 1; this helps us treat this term implicitly in our numerical methods and thus strictly conserve 
mass. 

In summary, concentration fluctuations are modeled using the equations, in addition to incompressibility, 

V;r=TjV2v + V. (y2t7k^ir) (5) 

dfC + v-Vc =D'V ■ (Vc + cSt^T) . (6) 


The boundary conditions on the top and bottom boundaries (sapphire plates) are zero flux boundary conditions, giving 
the Robin boundary condition Vc = —cSj^T at the boundaries. The initial condition we start from is a uniform solution 
of concentration co = 0.018; with time this decays to an exponential average profile that solves ^ (c) = — {c)St'^T (see 
Fig. 1). In our simulations we have accounted for the initial transient in establishing the concentration profile across 
the sample. Based on measurements of the time response of the PID servos that control the temperature of the sapphire 
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windows, the temperature gradient in the y direction is modeled with the following empirical fit as a function of time, 
Vr = ^ [l — exp (0.00540f — 0.0602?*'-^^®)] y, where AT = 13.25 is the estimated steady-state temperature difference. 

The spatial discretization are essentially identical to those in our previous work.^^ The simulations of the transient 
development of concentration fluctuations used the overdamped temporal integrator summarized in Algorithm 3 in Ref. 24. We 
perform fully three-dimensional simulations on a domain of dimensions 0.654 x 0.1 x 0.654 (this tries to match the smallest 
wavenumber in the simulations with the wave numbers measured with CCD camera in the experiments), discretized on a 
256 X 40 X 256 grid, using a time step size of At = 10 s. The structure factors S were averaged radially to obtain 

S {k-,t), where k = ^k^ + k^, using an averaging procedure that mimics that used in the analysis of the experimental data. Note 
that in this case it is possible to obtain the same results using two-dimensional simulations (k^ = 0) because of the symmetries 
of the linearized equations. Nevertheless, we chose to obtain three-dimensional results directly comparable to experiments. 
Sixteen independent simulations were performed and the results averaged to reduce statistical noise and estimate statistical 
error bars. To obtain the static structure factor at steady state, we initialized the system using the steady state concentration 
profile, and fixed VT = AT jL. For these steady-state runs we used a time step size Af = 80s and averaged over a single run of 
2000 time steps (corresponding to about 44 hours of physical time) skipping the initial 200 time steps in the analysis in order to 
allow the system time to reach a statistical steady state. 

Temperature fluctuations 

The dynamics of velocity and temperature (2,4) occur at similar time scales and must be integrated together; it is not justified to 
delete the inertial term pdtV in the velocity equation as it was for concentration. Therefore, for temperature we solve the system 
of equations 

pd,v + 'V7t ^ri'V^v + 'V- (^^JlriksToW^ (7) 

dj + vVT ^kV^T. (8) 

The boundary condition for temperature at the top and bottom walls are Dirichlet conditions, with T{y = Q,t) = 304.6 at one 
of the boundaries, and T{y = L,t) = 291.4 at the other wall, leading to a linear steady state temperature profile. Since for 
temperature we are not interested in the transient behavior, but rather only the steady state static structure factor, the initial 
temperature field is set to be the linear steady state. 

The spatial discretization is identical to that for concentration, in fact, our computer code does not distinguish between 
temperature and concentration since the equations are essentially identical. The temporal integrator is the inertial scheme 
summarized in Algorithm 1 in Ref. 24, requiring a much smaller time step size Af = 0.0016 s in order to resolve the fast 
vorticity dynamics. In this case we perform two dimensional simulations in a domain of dimensions 0.8 x 0.1 on a grid of 
256 X 32 grid cells. We average over 16 simulations of 5 • 10^ time steps each (corresponding to about 800 s of physical time), 
skipping the initial 5 10'* time steps. 
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